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Abstract 

We present a theoretical study of ghost imaging by using blackbody radiation source. A 
Gaussian thin lens equation for the ghost imaging, which depends on both paths, is derived. 
The dependences of the visibility and quality of the image on the transverse size and 
temperature of the blackbody are studied. The main differences between the ghost imaging 
by using the blackbody radiation and by using the entangled photon pairs are image-forming 
equation, and the visibility and quality of the image. 

PACS numbers: 42.30. Va, 42.25.Kb, 42.50.Ar, 42.25.Hz 

Ghost imaging were firstly reahzed by using entangled photon pairs generated in 

spontaneous parametric down conversion in 1995 [1,2]. The name ghost comes from the 

facts that an object in one path produces an image in another path in the measurement of 

coincident counting rates, and the image depends on both paths. Since then, many theoretical 

and experimental studies on this subject have been published [3-10]. Recently, there were 

discussions about whether quantum entanglement is necessary in the ghost imaging and 

whether a ghost imaging experiment can be realized with classical source [11-14]. Bennink et 

al. presented their classical ghost imaging and interference experiments [12-13]. However, 

they did not give the imaging formation equation. Furthermore, Angelo et al. [14] claimed 

that the classical ghost image [12] is a shot by shot, point-to-point projection. Gatti et al first 

pointed out theoretically that the ghost imaging can be achieved with truly incoherent hght 

[15], and then they worked under a particular path configuration, which is not the same as 

used in the quantum ghost image and interference experiments [1,2], and consequently they 

could not give the image formation equation and the interference fringe equation [15]. More 

recently, Cheng et al study the coincidence interference with a complete incoherent light by 

using classically theory [16], and as a matter of fact, the visibihty of the fringes in their case 

is zero, that is to say, there are no interference fringes at all. In this paper, we study the ghost 

imaging by using a blackbody radiation source [17] and derive the image formation equation, 

and investigate the visibihty and quality of the image and their dependence on the 
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temperature and the surface size of the blackbody. The difference between the imaging with 
the entangled photon pairs and the imaging with blackbody radiation is also discussed. 




FIG.l The scheme for ghost interference with blackbody radiation 

The scheme for the ghost imaging with blackbody radiation is shown in Fig.l. The 
blackbody radiation first propagates through a beam split, then propagates through path one 
and two to detector one and two (Dl and D2). In path one (from the blackbody to Dl), 
between the beam split and Dl, there is an object (a double-sIit aperture) with its transmision 
function H{v). In path two (from the blackbody to D2), there is a lens with focal length / 
between the beam split and detector two, and the distances from the blackbody to the lens 
and the lens to D2 are and I2, respectively. The scheme is the same as in [I] except that the 
lens in parth two in order to emphasize the object and the image are in different paths. (The 
lens in path one will be discussed at the end.) The coincident counting rate is proportional to 
the fourth order correlation function, g'^' (mj , ^2 ) . 

According the second order and fourth order optical coherence theory, the fourth order 
correlation function between the two detectors obeys the following integral formula for the 
blackbody radiation [17], 



2 



G''' (m, , m J = (e{m, )e{m2 )E* {u^ )E* (m, )) 
— 1 1 1 (''^1 ■> Wi )hl , Mj )/i2 , Mj )/i2 , Mj )x ^^(x, )£'(jf2 (^3 )E* {x^ ^dx^dx2dx.^dx^ 



= </(Ml)></(M2)> + |r(Mi,M2)| , (1) 

where h^{x^,u^), /i2(^2'"2) the response functions of the two paths through which the 
blackbody radiation passes, and 



/(M; ) >= J J^i (^1 5 M,- (^2 ' )-^' (^2 fjdx^dx2 i = 1,2 , 



(2) 



r(Mj,M2)- J° (£'(jCi)£'*(jC2))^lUl'"l)^2(^2'"2)'^l'^-^2 • (3) 

< /(m, ) > is the second order correlation function at the same space point (the intensity at the 
i-th detector ), and depend only on the i-th path. r(M^,M2) is the second order cross 
correlation function at two different points, which is related to both detectors and depends on 
both paths. Both</(Mj)> <I(u2)> and r(Mj,M2) have contribution to the coincident 
counting rates. In Eq. 1, the fourth order correlation is replaced by the second order 
correlation under the condition of< E(x.) >= for the blackbody radiation [17]. 

For the blackbody radiation the second order correlation function takes the following 
form [17]: 
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exp[«k • (r2 - ^k. 



(4) 



le^ )^ t^^{h(ol kgT)-\ 

where T is the blackbody' s temperature, and kg is the Boltzmann's constant. For simplicity, 
we consider the source in one-dimension. The numerical calculation shows that 
{e{x^E*{x2^ is a quasi-Gaussian distribution over (jCj-jc2) with a temperature dependent 
width [17]. Therefore, we can approximately write the second order correlation function as 



( +x^^ 
£'{xi)£'*{x2))~/(Xi,C7 )exf ' ^ 



exp 



(X[ X2 ) 



(5) 



where cr^ is the correlation length and is temperature dependent. As the temperature goes to 
infinity, (7 approaches zero and I{r^,<J ) infinity. Here we introduced a Gaussian function 
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exf 



( 2 2\ 



to take into account the finite surface size of the blackbody with <7^ the 



transverse size. 

With the help of Collin's formula and the detail information of the two paths, h^{x^,u^) 
and /i2(^2'"2) [18], we obtain 

</(w,)>=V— £ r r [<E(x,)E(x,)>H{v,)H'{v,)exT^ -2x,v,+v^) 

A Z^Z^ ^ J-oo J-oo J-oo 2v 



i7t I 2 ^ i\ 
xexp — - — -zViMj +Mi j 

L ^2 

6) 



^^^^[X2-2X2V2+V2) 
Az, 



Az^ 



exp -^{yl -2V2M1 +uf) dx^dx2dv^dv2 
Az-, 



For a fixed Mj, < /(Mj) > is a constant, and it is easy to show that < /(Mj) >= fio=^<^^^t. The 
image pattern is determined by the cross correlation function, 

r(Mi , M2 ) =< e{u^ )e* ) > 

r ■ 

r r r < E(x, )E* {x^ ) >H{v, ) exp - — (xf - 2jc, V, + vf ) - — (vj' - 2v,u, + ) 
j-00 j-00 ^ 



exp ^^{ujxl - 2X2U2 +d2ul) 



Ab. 



exp[- +22) + '^(^1 ~ h )]dXidx2dVi 



(7) 



where fl,,&2'<^2^'^d rfj^re the optical transfer matrix elements between the blackbody and 
detector two. 
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(8) 



If we set 



Eq. (8) is reduced to 



1 1 1 

- + — = 



K - Zi I2 f 



(9) 



«2 K 
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(10) 
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By substituting Eqs. (5) and (10) into Eq. (7), we have the cross correlation function. In order 
to obtain an analytical expression, we assume blackbody is sufficiently large ( (7^ — > oo ) and 
at very high temperature (a — > 0), and we have 



j2 



H 



f \ 

V^2 J 



(11) 



This is a perfect image of the object with an amplification of , which is obtained under Eq. 
(9), the Gaussian thin lens equation. We emphasize that the conditions ( (7^ — > and 
a J — > oo) are only for obtaining the analytical expression of Eq. (11) . When Cj andcr^ are 

finite or Eq. (9) is not satisfied, the image needs to be obtained numerically. In Fig. 2, we plot 

the evolution from the ghost image into an interference pattern of a double slits with slit 

width , and distance of the two shts d, when we vary from satisfying to not satisfying Eq. 

(9). The transmission function of the double slits is 

da d a , d a d a ^ ^ . 

H(v) = 1 tor <v< h — and <v< — h— and =0 otherwise, 
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FIG.2 The image pattern of a double sht in the scheme of Fig.2 with blackbody radiation 
for different l^- (a) = 20mm (imaging case), (b) = 20.5mm, (c) = 21.5mm with 
A = 102nm , a = 0.01mm , d = 0.03mm , Zj = 10mm , Z2 = 40mm , f = \Omm , 
/j = 30mm and <7^ = 5mm , a - 0.00001 . 



Now we need to consider the visibihty of the image. We define the visibility of 
the image as 
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Max\r{u^ — OjMjl 



(12) 



Max{l{u, =0))(/(m2))' 



Under the condition of obtaining Eq. (11), — > and cr^ — >oo, the visibility of the 

image is zero. The images of the double slit aperture for different transverse sizes and 
transverse coherence widths (temperature) of the blackbody are shown in Fig.3 and Fig. 4, 
respectively. From Fig.3, we can find that when the source's transverse size cr^ increases, 
the quality of the image increases, while the visibility decreases. From Fig.4, we can find 
that when the source's transverse coherence width <7^ decreases, quality of the image 

increases, while the visibihty decreases. 



Fig.3 The image of the double slit aperture for different source's transverse sizecr^ : (a) 
0.1mm, (b) 1mm and (c) 5mm with A = 102nm,a - 0.01mm, d - 0.03mm , - 10mm , 
Z2 = 40mm , / = 10mm , = 30mm and /j = 20mm , a =0.001mm. 
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Fig.4 The image of a double slit aperture for different source's transverse coherence 
width £7^ : (a) 0.00001mm, (b) 0.0005mm and (c) 0.003mm with A = IQlnm , a = O.Olmm , 

d - 0.03mm , = 10mm , - "^Omm , f = 10mm , = 30mm and = 20mm , 
<Tj =5mm. 



Now we study the dependence of the visibility and quality of the image on the source's 
transverse coherence width (T^ and transverse sizecr^ in details. We define an image quality 
factor, 

|r(M, = 0,1/2)1^ 

r(Mj = 0,M2) 



-\H{u2)\^du 



(13) 



r\H{u^)\^ du^ 

J-00 

Small Q value corresponds to high image quality. The dependences of the visibility and 
quality of the image of the double slit aperture on the transverse coherence width are shown 
in Fig. 5 and Fig. 6, respectively. High quality is companied by poor visibility, and good 
visibility companied by low quality. In order to observe the classical ghost image with 
blackbody radiation, the selection of suitable transverse size and transverse coherence width 
is essential. 
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Fig.5. Evolution of the visibility of the image of a double slit aperture versus transverse 
coherence width cr^ for different transverse sizes of the source, (7j . (a) Imm (b) 5mni (c) 

10mm with X = IQlnm , a = 0.01mm , d = 0.03mm , Zj = 10mm , Z2 = 40mm , / = 10mm , 
/j = 30mm and 1 2 - 20mm . 
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Fig.6. Evolution of the quahty of the image of a double slit aperture versus transverse 
coherence width cr^ for different transverse sizes of the source, <Jj . (a) 1mm (b) 5mm with 

A = 702nm , a = 0.01mm , d = 0.03mm , = 10mm , Z2 = 40mm , / = 10mm , /j = 30mm 
and I2 - 20mm . 
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Now, we ask ourselves whether the classical ghost unage is the same as the quantum 
ghost image with the entangled photon pairs. Comparing Eqs. (9) with the corresponding 
equations for the quantum ghost image [2], the following differences have been found: (1) 

/j + Zi in the imaging formation equation for the quantum case is replaced by /, - Zj for the 
blackbody. (2) For li> z^, we have an inverted image and enlarged image with > 2/ or an 
inverted and reduced image with 2/ > > / ; for li< z^, we have an erect and reduced image 
with < /2 < / • For the quantum ghost image, the image is always inverted. (3) Another 
important difference between quantum and classical ghost image is the visibility. In the 
Quantum ghost image, good quality and high visibility can be achieved simultaneously [1, 2], 
but this is impossible in the ghost image with blackbody radiation. The differences come from 
the classical and quantum nature of the correlation functions. In quantum case, the cross 

correlation function is r(jCj,JC2)°« {|0,0|£'(xj )£'(jC2 )|^^ due to the entanglement, while in the 

classical case (no entanglement) it is Eq. (3) where we have E*(x2) instead of ^'(jCj). If a 
50% phase conjugate mirror is used for the beam spht, we will have Qi - Zq) + (z^ - z^) 
instead of (l^ - Zq) - (Zi - Zq) with Zq the distance between the black body and the beam spht. 

If the lens is in path one, we will have — h — = — with 5, and 5, the distances from the 

S,-S, S, f 

blackbody to the lens and to D2, and ^3 from the lens to the object. (In quantimi case, - ^2 is 

replaced by + ^2 .) Further more, in quantum case we have G^^\x^,X2) =\ r(x|,X2 )l' , that IS 
to say , there is no background noise and consequently the visibility is high. In the classical 
case the background noise limits the visibility. 

If we take out the lens in path two (see Fig. 1), we have the ghost interference. With the 
same calculation above, we can obtain the same interference fringes as obtained from the 
quantum ghost interference except the replacement of - Zj by + Zj . 

In conclusion, we have invested the ghost image created with blackbody radiation by 
using the optical coherence theory. The ghost image formation depends on both paths. To 
obtain the ghost image, a Gaussian thin lens equation must be satisfied. The ghost image is 
gradually blurred out, when the temperature (or the size) decreases. The quality of the ghost 
image increases with the increase of the temperature and the size of the blackbody. As the 
temperature and the size both increase to infinity, we will have a perfect image, but the 
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visibility goes to zero. The nature of the ghost imaging is due to the entanglement in the 
entangled photon pair and is due to the Hanbury Brown-Twiss Effect (low coherence with 
fluctuations but not completely no coherence) in the blackbody radiation. 
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